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Abstract — In this paper we study the capacity region of the 
multi-pair bidirectional (or two-way) wireless relay network, 
in which a relay node facilitates the communication between 
multiple pairs of users. This network is a generalization of the 
well known bidirectional relay channel, where we have only one 
pair of users. We examine this problem in the context of the 
deterministic channel interaction model, which eliminates the 
channel noise and allows us to focus on the interaction between 
signals. We characterize the capacity region of this network when 
the relay is operating at either full-duplex mode or half-duplex 
mode (with non adaptive listen-transmit scheduling). In both 
cases we show that the cut-set upper bound is tight and, quite 
interestingly, the capacity region is achieved by a simple equation- 
forwarding strategy. 

I. Introduction 

Cooperative communication and relaying is one of the 
important research topics in wireless network information 
theory. The basic model to study this problem is the 3-node 
relay channel which was first introduced in 1971 by van der 
Meulen [1] and the most general strategies for this network 
were developed by Cover and El Gamal [2]. While the main 
focus so far is on the one-way-relay channel, bidirectional 
communication has also attracted attention. Bidirectional or 
two-way communication between two nodes was first studied 
by Shannon himself in [3]. Nowadays the two-way communi- 
cation where an additional node acting as a relay is supporting 
the exchange of information between the two nodes (or one 
pair) is gaining increased attention. Some relaying strategies 
for this one-pair two-way relay channel, such as decode- 
and-forward, compress-and-forward and amplify-and-forward, 
have been analyzed in [4]. Network coding type techniques 
have been proposed by [5], [6], [7], [8] (and others) in order 
to improve the transmission rate. 

The two way relay channel problem can be generalized 
to a multi-pair (or multiuser) setting in which the relay 
facilitates the communication between multiple pairs of users. 
In [9] authors analyzed the case that the relay orthogonalizes 
different two-way transmissions by a distributed zero forcing 
algorithm and then multiple pairs communicate with each 
other via several orthogonalize-and-forward relay terminals. 
In [10], [11] authors investigated this problem for interference 
limited systems in which each pair of users share a common 
spreading signature to distinguish themselves from the other 
pairs, and proposed a jointly demodulate-and-XOR forward 
strategy. However, so far no attempt has been done to char- 



acterize the capacity region of this network, and the optimal 
relaying strategy is unknown. 

In this paper we study the information theoretic capacity 
of the multi-pair bidirectional wireless relay network. We 
examine this problem in the context of the deterministic 
channel interaction model. The deterministic model, studied 
by Avestimehr, Diggavi, and Tse [12], simplies the wireless 
network interaction model by eliminating the noise and allows 
us to focus on the interaction between signals. This approach 
was successfully applied to the relay network in [12], and 
resulted in insight in terms of transmission techniques which 
further led to an approximate characterization of the noisy 
wireless relay network problem [13]. This approach has also 
been recently applied to the bidirectional relay channel prob- 
lem [14], which again resulted in finding near optimal relaying 
strategies as well as approximating the capacity region of the 
noisy (Gaussian) bidirectional relay channel. 

Inspired by these results, in this paper we apply the deter- 
ministic model to the multi-pair bidirectional relay network, 
and analyze its capacity when the relay is operating at either 
full-duplex mode or half-duplex mode (with non adaptive 
listen-transmit scheduling). In both cases we exactly character- 
ize the capacity region and show that the cut-set upper bound 
is tight. Quite surprisingly, we show that the capacity region is 
achieved by a simple equation-forwarding strategy, in which 
different pairs are orthogonalized on the signal level space and 
the relay just re-orders the received equations created from 
the superposition of the transmitted signals on the wireless 
medium and forwards them. 

The paper is organized as follows. In Section HT1 we state the 
precise definition of the problem. We will present our main 
results in Section|III]and characterize the exact capacity region 
of the full-duplex and half-duplex multi-pair deterministic 
bidirectional relay networks. Then in Section [IV] we illustrate 
further structures in the optimal relaying strategy and finally 
conclude in Section [V] 

II. System model 

The system model for the M-pair bidirectional relay 
network is shown in Figure Q] In this system M pairs 
(Ai, Bi), . . . , (Am,Bm) aim to use the relay to communicate 
with each other (i.e. A\ and B\ want to communicate with 
each other, and so on). The relay can operate on either full- 
duplex or half-duplex mode. In the full-duplex mode it is able 




(a) Uplink 




Fig. 1. 
network. 



(b) Downlink 

The system model for M pair bidirectional deterministic relay 



to listen and transmit at the same time, while in the half- 
duplex mode it can only listen or transmit at a particular time. 
In the half-duplex scenario, we only consider the case that the 
listen-transmit scheduling is non-adaptive and the relay listens 
a fixed t fraction of the time and transmits the rest. Although 
t can not change adaptively as a function of the channel gains, 
one can optimize over t beforehand. 

We use the deterministic channel model to model the 
interaction between the transmitted signals. The deterministic 
channel model was introduced in [12]. Here is a formal 
definition of this channel model. 

Definition 2.1: (Definition of the deterministic model) 
Consider a wireless network as a set of nodes V, where \V\ = 
N . Communication from node i to node j has a non-negative 



integer gai 



associated with it. This number models the 



channel gain in a corresponding Gaussian setting. At each time 
t, node i transmits a vector x,[t] G F| and receives a vector 
Yi[t] S F2 where q = maxj^n^ 3 -)). The received signal at 
each node is a deterministic function of the transmitted signals 
at the other nodes, with the following input-output relation: if 
the nodes in the network transmit xi[i],x 2 [i], . . .xjv[i] then 
the received signal at node j, 1 < j < N is: 

N 

y^] = ^S«-"«x fe [i] (1) 

k=l 

for all 1 < k < N, where S is the q x q shift matrix and the 
summation and multiplication is in F2. 

Now that we have defined the deterministic channel model 
we can apply it to the multi-pair bidirectional relay network. 
A pictorial representation of an example of such network with 
two pairs is shown in Figure |2] In this figure each little circle 
represents a signal level and what is sent on it is a bit. The 
transmit and received signal levels are sorted from MSB to 
LSB from top to bottom. The channel gain between two nodes 
i and j indicates how many of the first MSB transmitted signal 
levels of node i are received at destination node j. Now as 

'Some channels may have zero gain. 
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(a) Uplink 



(b) Downlink 



Fig. 2. The pictorial representation of a two-pair bidirectional deterministic 
relay network with channel gains n^jj = 3, ng^ = 2, «a 2 -R = 2, 
n-B 2 R = 1. n RAl = 2, n RBl = 3, n RA2 = 1 and n RB2 = 2. 



described in the channel model (Q]i, at each received signal 
level, the receiver gets only the modulo two summation of the 
incoming bits. 

III. Capacity region of the multi-pair 

BIDIRECTIONAL DETERMINISTIC RELAY NETWORK 

In this section we study the capacity region of the multi- 
pair bidirectional deterministic relay network. First we state 
our main result and then in the rest of this section prove it. 

Theorem 1: The capacity region of the multi-pair bidi- 
rectional deterministic relay network both (full-duplex and 
half-duplex), described in Section [TTJ is equal to the cut-set 
upper bound region, and it is achieved by a simple equation 
forwarding relaying strategy. 

For example in the case that we have only two pairs and the 
relay is operating on the full-duplex mode, the cut-set upper 
bound on the capacity region is given by, 



Ra 1 b 1 <mm(n AlR ,n RBl ) 
Rb 1 A 1 < min (n BlR ,n RAl ) 
Ra 2 b 2 < min (n A2 R,n R B 2 ) 
Rb 2 a 2 < min (n B2R ,n RA2 ) 



(2) 
(3) 
(4) 
(5) 



R AlBl + Ra 2 b 2 < min (max (n AlR ,n A . 2R ) ,max (n RBl , n RB2 )) 



(6) 



Rb 1 a 1 + Rb 2 a 2 < min (max (n Bl R, t»b 2 k) , max (n RAl , n RA2 )) 



(7) 



Ra 1 b 1 + Rb 2 A 2 < min (max (n Al R, "B 2 «) > max (n RBl ,n RA2 )) 



<-B 1 A 1 



(8) 

Ra 2 b 2 < min(max(n Blfl ,n A2fl ) , max (n RAl , n RB2 )) 

(9) 

Now consider the network shown in Figure [2] It is easy to 
check that the rate tuple (R AlBl , Rb^, Ra 2 b 2 , Rb 2 a 2 ) = 
(2, 1,1,1) is inside its cut-set region. In Figure[3]we illustrate a 
simple scheme that achieves this rate point. With this strategy, 
the nodes in the uplink transmit 



XAi 



[01,1,01,2,0]*, x Bl = [61,1, 0,0]* 



xa 2 = [0,a 2 ,i,0] , x Bl = [62,1,0,0] 
and the relay receives 

2/R=[oi,i, ai,2©&i,i, a 2i i©6 2i i] t 



(10) 



Then the relay will re-order the received equations and trans- 
mit 

xr = [02,1 © &2,i, ai, 2 ©&i,i, ai,i]' (11) 




(b) Downlink 



Fig. 3. The scheme that achieves rate point (2, 1, 1, 1). 

Then node A\ receives 1x1,2 © 61,1 and since it knows 01,2 
can decode similarly node B\ can decode ciij and a\ t %, 
node A2 can decode 62.1 and finally node B2 can decode 02.1. 
Therefore we achieve the rate point (2, 1, 1, 1). 

There are some interesting points about this particular 
achievability strategy, 

• There is no coding over time, 

• There is no interference between different pairs on the 
same received signal level at the relay, 

• The relay just re-orders the received equations and for- 
wards them. 

We call a scheme with these properties an equation forwarding 
scheme. A natural question is whether we can always achieve 
any rate point in the cut-set bound region with an equation- 
forwarding scheme? 

Quite interestingly, in the next lemma we show that indeed 
this is always possible. 

Lemma 1: Any integral rate tuple inside the cut-set upper 
bound on the capacity region of the full-duplex M-pair bidi- 
rectional deterministic relay network is achieved by a simple 
equation forwarding scheme. 

Proof: Due to lack of space, we only prove the lemma 
here for M = 2 pairs, but the proof can be easily gen- 
eralized to any arbitrary number of pairs. We use induc- 
tion on the sum rate to show that every integer 4-tuple 
(Ra 1 b 1 ,Rb 1 a 1 , Ra 2 b 2 ,Rb 2 A 2 ) satisfying the cut set bound 
is achievable by equation forwarding. For convenience we 
consider two separate cases, and mention how in each case 
we can assign a relay signal level in uplink and downlink for 
serving a session worth 1 bit. We then show that the reduced 
rate tuple is within the cut-set region of the reduced network. 

2 with integer components. 



We use an example along the proof to illustrate the ideas. The 
example network is shown if Figure [5] and one can check that 
the rate-tuple (3,1,2,1) is in its cut-set region. 

Case 1: There is a pair where both nodes have nonzero 
transmission rates. Without loss of generality we may assume 
that Ra 1 b 1 and RbiAi are both nonzero. Our goal is to assign 
one up-link signal level and one down-link signal level to the 
(Ai,Bi) session. A\ and B\ will then transmit one bit at the 
assigned uplink level to the relay, and the relay will transmit 
the received equation at the assigned down-link level to both 
A\ and B\. After doing so and removing the assigned levels, 
the network will be converted to an equivalent network with 
lower capacities and reduced rates. If we next show that the 
new rates satisfy the cut-set bounds of ((2j to © for the new 
network, then by induction we can say that the cut-set bound 
is achievable in this case. 

We claim that in up-link the highest signal level connected 
to both A\ and By, and in down-link the lowest signal level 
connected to both A\ and B\ are the appropriate levels that 
we are looking for. We refer to these levels by l u and Id 
respectively. For the sample network of Figure [5] l u and Id 
correspond to the third relay's circle from the bottom in UL, 
and third relay's circle from the top in DL respectively. Now 
the claim is after removing l u and Id from UL and DL, the 
new rate tuple {Ra 1 b 1 — 1, Rb^ — L Ra 2 b 2 , Rb 2 a 2 ) satisfies 
the cut-set bounds of the new network. First, note that by 
removing the levels l u and Id, each of the values Ua^r, nBiR, 
riRAi and urb 1 will decrease by exactly 1, and the other 
four capacities Ua 2 r, ub 2 r, Ura 2 and urb 2 will decrease 
by at most 1, depending on the situation. So, after removing 
the levels and reducing Ra 1 b 1 and Rb 1 a 1 by 1, equations 
© and © clearly continue to hold . Besides, equations © 
to © all remain valid, since in each of them the L.H.S is 
decremented by 1 whereas the R.H.S by at most 1. Because 
of the symmetry, we only have to resolve the case of equation 
©. Let's say © is violated because of eliminating the uplink 
signal level l u . This means prior to this we should have 
had: Ra 2 b 2 < nA 2 R and on the other hand l u is connected 
to A2. Because of the choice of l u , the latter means that 
riA 2 R > min{n Al R,n Bl R}. Let's say tla 2 r > n Al R. This 
and © imply that iia 2 r > Ra 1 b 1 + Ra 2 b 2 - But we already 
have Ra 2 b 2 = n A 2 R, which means Ra 1 b 1 = 0, which 
contradicts our initial assumption. Similarly if ua 2 r > ns 1 R, 
from © we conclude Ua 2 r > Rb-lAi + Ra 2 b 2 - This and 
the fact that Ra 2 b 2 = n A 2 R, imply Rb ± Ai = 0, which again 
cannot be true because of the original assumption. 

In a similar way we can show that the removal of Id can't 
spoil validity of ©. The reason is if Ra 2 b 2 = tirb 2 an d 
Id is in the transmission range of relay to B2 (i.e. the circle 
corresponding to level Id in Figure [5] is connected to B2) then 
riRB 2 > vam{nRAi,nR Bl }. Now if n RB2 > n RAl then from 
© nRB 2 > Ra 2 b 2 + Rb 1 a 1 - But Ra 2 b 2 = n R B 2 which 
means RbxAx = 0. Likewise, if jirb 2 > fiRBi from equation 
© we conclude urb 2 > Ra 1 b 1 + Ra 2 b 2 , which implies 
Ra 1 b 1 = 0, which is a contradiction. 



Let's apply this scheme to the example network. We should 
first take the (^4i,i?i) pair and serve them through one signal 
level in UL and one level in DL. Next, we take the (A 2 ,B 2 ) 
and similarly assign corresponding levels in UL and DL to 
them. These two steps are shown in Figures and |4j). For 
clarity, the removed signal levels are dotted in each step. The 
remaining unserved rates are (2,0,1,0). 

Case 2: Every session has a node with zero rate. W.L.G., 
say Rb 1 a 1 — Rb 2 a 2 = (only A\ and A 2 want to 
communicate with their end parties). We show that A\ and A 2 
can respectively transmit Ra 1 b 1 and Ra 2 b 2 bits at different 
signal levels to the relay, and the relay can forward them to 
B\ and B 2 by putting on distinct signal levels in down-link 
without any interference. In fact, if ua x r < nA 2 R, then A\ can 
put Ra 1 b 1 bits on the lowest signal levels in uplink and A 2 
can put Ra 2 b 2 bits on the levels on top of it. This is feasible 
because of equations (fJJ and ©. However, if Ua 2 r < Ha^Ri 
A 2 can put his bits on the lowest levels and A\ will transmit on 
top of it. Similarly, in DL if urbi < ^/?b 2 the bits addressed 
to B\ will be put on the highest signal levels and the bits 
addressed to B 2 on the signal levels below those of B\ . Figure 
@fc shows how this idea is applied to our example network. 
The final equation forwarding configuration that achieves the 
rate-tuple for this example is shown in Figure [5] ■ 

Now to prove the main Theorem Q] for the full-duplex 
scenario, we just need to show that all corner points of the 
cut-set bound region are achieved by an equation forwarding 
scheme. Note that since all coefficients of the hyperplanes of 
the cut-set bound region are integers, then all corner points 
of the region must be fractional. If a corner point R is 
integral then by Lemma Q] we know that it is achieved by an 
equation forwarding scheme. If it is not integral then choose 
a large enough integer Q such that QR is integral. Now note 
that Q instances of a deterministic network overtime, is the 
same as the original network except all channel gains are 
multiplied by Q. Now since QR is integral and is obviously 
inside the cut-set upper bound of the big network (where all 
channel gains are multiplied by Q), then by Lemma Q] it is 
achievable by an equation-forwarding scheme. This scheme 
can be simply translated to an equation forwarding scheme on 
the original network over Q time-steps. Therefore the corner 
point Qj- = R is achievable. 

Similarly we can also prove Theorem Q] for the half-duplex 
scenario. In this case relay listens t fraction of the time 
and transmits the rest. Without loss of generality assume 
t is a fractional number (otherwise consider the sequence 
of fractional numbers approaching it). Then choose a large 
enough integer Q such that Qt is integer. Then consider Q 
instances of the network over time, such that for Qt instances 
the relay is listening and in the other Q(l — t) instances it is 
transmitting. After concatenating these instances together, the 
resulting network can be thought of as a full-duplex multi- 
pair network where the uplink channel gains are multiplied 
by Qt and the downlink channel gains are multiplied by 
(1 — t)Q. It is easy to verify that the cut-set bound region 



of this network is just the cut-set bound region of the original 
half-duplex network expanded by Q. Now by Lemma Q] and 
the previous argument, we know that the capacity region of 
this full-duplex multi-pair bidirectional network is equal to its 
cut-set upper bound region and is achieved by an equation 
forwarding scheme. Now note that any equation forwarding 
scheme in this full-duplex network can be translated to an 
equation forwarding scheme in Q instances of the original 
half-duplex network; Qt instances the relay is in the listen 
mode to get the equations and (l—t)Q instances in the transmit 
mode to send the equations. Therefore the cut-set upper bound 
region is achievable and the proof is complete. 

IV. Remark 

Although the scheme we provided is an inductive way of 
level assignment and seems unstructured in the sense that it 
assigns signal levels on a greedy basis, one can say more 
about these assignments using certain observations. First of 
all, note that in this equation forwarding scheme we have in 
general 2M types of equations that the relay might decode. 
Namely, M types of equations getting bits from one user 
of a session, and M types of equations getting bits from 
both users of the same pair. Refer to the example network 
of Figures @] and [5] and observe that in the final config- 
uration all equations of the same type were concatenated 
together both in UL and DL. In general one can serve all 
equations of the same type at once by choosing a pair with 
nonzero rates and serving them until one of the rates is 
zero. For example, assuming RaxBi < Rb-iAi an d Ra 2 B! < 
Rb 2 a 2 , instead of reducing (R Ai b 1 ,Rb 1 a 1 ,Ra 2 b 2 Mb 2 a 2 ) to 
{RAtBt - - 1,Ra 2 b 2 ,Rb 2 a 2 ) one can reduce it to 

(0,Rb 1 a 1 — Ra 1 b 1 ,Ra 2 b 2 ,Rb 2 a 2 ) all at once. Then the same 
thing can be done for the other pair. Now by rearranging the 
signal levels, it is easy to show that in the final configuration, 
all the equations of the same type are in concatenation. The 
conclusion of the above argument is that for any network con- 
figuration, one can always find 2M groups of disjoint signal 
levels (both in UL and DL) with the following properties: 



« Group 2i — 1: min{i? J 4 i B i , .Rb^} levels connected to 

both Ai and Bi. 
. Group 2i: max{i?^ iBi , -RemJ - min{i?A I s I , Rb^} 

levels connected to the one of Ai and Bi with higher 

transmission rate. 



Furthermore, levels of each group are all in concatenation. 
Intuitively this observation suggests that in the noisy case each 
user need to only break the transmit signal to at most 2M 
different levels, independent of the channel gain strengths. 

V. Conclusions 

In this paper we studied the multi-pair bidirectional relay 
network which is a generalization of the bidirectional relay 
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Fig. 4. Illustration of the inductive algorithm introduced in Lemma [T] 




(b) Downlink 



Fig. 5. Illustration of the resulting equation forwarding scheme of the 
inductive algorithm. 

channel. We examined this problem in the context of the deter- 
ministic channel model introduced in [12] and characterized its 
capacity region completely in both full-duplex and half-duplex 
cases. We also showed that the capacity can be achieved by 
a simple equation-forwarding strategy and illustrated some 
structures on the signal levels that these equations are created 
at. In ongoing work, we have made progress on using these 
insights to find an approximate capacity characterization of 



the noisy (Gaussian) version of this problem. We hope to 
completely answer this question in a future work. 
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